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Abstract 

We present an explicit convergent iterative solution for the lowest energy state 
of the Schroedinger equation with a generalized double well potential V = \(x 2 — 
l) 2 (x 2 + a). The condition for the convergence of the iteration procedure and the 
dependence of the shape of the groundstate wave function on the parameter a are 
discussed. 
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1. Introduction 



This paper is stimulated by an interesting question raised by Roman Jackiw[l] 
concerning the extent of validity of the convergent iterative method[2-5] that 
we have developed for the iV-dimensional generalization of the double well 
potential 

V = |V - rlf. (1.1) 

More specifically, whether the method is equally applicable to a different 
sombrero-shaped potential 

V = ^(r"--rlf^ + 2rl). (1.2) 
The latter has several new features[l,6]. When g 2 = 1 and 

4 = 1(2 + N), (1.3) 

the groundstate wave function ip is simply 

e^ 4 (1.4) 

which peaks at r = 0. Yet for g 2 » 1, the maximum of ip has to be near 
r = ro. In particular, for a one-dimensional problem, as g 2 increases from 1, 
there would be a critical point when ip changes from having a single maximum 
at the origin to one with double peaks. Thus, it is of interest to examine 
whether our convergent iterative method works for g 2 = 1 as well as for g 2 
larger than 1. The purpose of this paper is to show that this is indeed the 
case. 

To simplify our discussions, we examine only the one-dimensional case in 
this paper. Let ip(x) be the groundstate wave function of the Schroedinger 
equation 

-~ + V(x)]tP(x)=E7P(x) (1.5) 
with 

V ( x ) = g -(x 2 -l) 2 {x 2 + a) (1.6) 
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and 



a > 0. 



(1.7) 



We note that in one dimension, the potential (1.2)-(1.3) is a special case of 
(1.6) with a = 2. For convenience of nomenclature, in dimension > 1 we call 
(1.1) the sombrero potential, and (1.2) the generalized sombrero potential; in 
one dimension, we call 



the double well potential and (1.6) the generalized double well potential. 

In Section 2, we give a brief review of our convergent iterative method for 
the potential (1.6). To ensure rapid convergence, we have established a rather 
effective theorem, called Hierarchy Theorem, provided a certain inequality can 
be satisfied. This inequality is proved in Section 3 for a = 2 (i.e., the potential 
(1.2) in one-dimension), and in the Appendix for an arbitrary positive a > a c = 
0.664. Some of the pertinent numerical results are given in Section 4. 

2. Trial Function and Iterative Equations 

To construct a good trial function (j){x) for the groundstate wave function, we 
follow the same steps developed for the double-well potential[2-5]. The main 
components of 4>{x) are constructed by using the Schroedinger equation (1.5) 
and extracting the first two terms of its perturbation series expansion of ip(x) 
and E in g^ 1 . 



2 v ' 



(1.8) 



Write 



(2.1) 



and set 



gS(x) = gS {x) + S x (x) + g 1 S 2 {x) + • • • . 

Correspondingly, 

E = gE + E 1 +g- 1 E 2 + ---. 
Substituting (2.1)-(2.3) into (1.5), we find 



(2.2) 



(2.3) 



S ! Q = (x 2 - l)Vx*Ta, 



(2.4) 
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and 



S[S' Q = ±SS - Eo. (2.5) 



Since the left side of (2.5) is zero at x = 1, 



E = isJKi) = VTT^. (2.6) 



Thus, 



, _ x{3x 2 + 2a - 1) - 2V1 + aVx 2 + a 

1 ~ 2(x 2 - l)(x 2 + a) ' ( j 
(Throughout the paper, prime denotes d/dx.) Next, introduce 

0+(ic) = exp[— gSo(x) — Si(x)} (2-8) 

and 

0_(x) = e^[-^ (-x) - 5i(a;)]. (2.9) 

We note that (f>+(x) and 4>-(x) are completely defined, except for a common 

arbitrary normalization factor. The trial function <j){x) is an even function of 
x, defined by 



4>(x) = 4>(—x) = 

with 



0+(x) + T 0_(x) for < x < 1, 

{1 + [T0_(l)/0 + (l)]}0 + (x) for x> 1, 



(2.10) 



r= 0+(O) = ga- Vl + a , . 

(f)'_(0) ga + VTT7i ' [ ' ' 



In the following, we assume 



1 



g > -VTT^ , (2.12) 



and therefore 



T > 0. (2.13) 

By construction, 0'(O) = 0. The functions 0(x) and 0'(x) are continuous 
everywhere. Since ip(x) and 0(x) are both even in x, we need only consider 

x > (2.14) 



in the following. 

By differentiation, <j) + and <j) satisfy two Schroedinger equations 

1 
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--<l>l + (V + u) ( J )+ = gE <j )+ 



and 

where 

and 
with 

Thus, for 



-\(j>" + {V + w)(j) = gE Q( j), 
u(x) = \(S[ 2 - S'l) 



w(x) = u(x) + g(x) 

q( x ) = J gE °i^kl for < x < 1 

x for x > 1. 



g'(x) = gE ( 



< x < 1, 



2.15) 

2.16) 
2.17) 
2.18) 

2.19) 

2.20) 
2.21) 

2.22) 

2.23) 
and therefore 

g'(x) < 0. (2.24) 

Thus, g(x) is a positive decreasing function of x. At x = 1, g has a discounti- 
nuity, decreasing to zero for x > 1. 

For the one-dimensional case of the generalized sombrero shaped potential 
(1.2)-(1.3), we have 

AT = 1, ro = l and a = 2. (2.25) 



From (2.8)-(2.9), it follows that 

t = e -2g5 (o) . e 2^ (x)_ 

0+ 

Within the range (2.20), 

'ty = (t) 2s(ji _ 1) v?Ti<o 



The corresponding potential is 

V(x) = ^(x 2 -l) 2 (x 2 + 2). (2.26) 

In the next section, it will be proved that for this potential, at all finite x > 0, 
we have 

u(x) > (2.27) 



and 

Thus 

and 



u\x) < 0. (2.28) 
w(x) > (2.29) 



w'(x) < 0. (2.30) 

Furthermore, w(oo) = 0. [The extension to the potential (1.6) with a^2 will 
be discussed in the Appendix.] 

Once the condition w'{x) < is established, we can apply the Hierarchy 
Theorem[3,4], as we shall discuss. 

Define 

H = -^^ + V(x)+w(x) (2.31) 

and write (2.16) as 

{H - gE )<f){x) = 0. (2.32) 

The Schroedinger equation we would like to solve is (1.5), which can be written 
as 

(H - gE )if>(x) = (w(x) - S)ij(x) (2.33) 

where 

6 = gE - E. (3.34) 
Multiplying (2.33) by 4> and (2.32) by ip. Their difference gives 

- = (w - E)U- (2-35) 

Hence, the ratio 

/ = (2-36) 



satisfies 

-\{4> 2 f)' = (w- £)4> 2 f- (2.37) 
Its integral over all x gives 

oo oo 

S = J W( j> 2 f dxj J (j) 2 f dx. (2.38) 



Eqs. (2.37) and (2.38) can then be solved by considering the iterative series 
{f n (x)} and {S n }, with 

-\(^a=(w-S n )^f n ^ (2.39) 

and therefore 

oo oo 

Sn = J W<t> 2 f n -idx/ J tffn-ldx. (2.40) 


We differentiate two different sets of boundary conditions: 

(I) / n (oo) = 1 for all n (2.41) 

or 

(II) / n (0) = 1 for all n. (2.42) 

Thus, in case (I) 

°r dv °r 

f n (x) = 1-2/ / (w(z) - S n )4> 2 {z)f n ^{z)dz (2.43) 
and correspondingly in case (II) 



v 



dy y 



f n (x) = 1 - 2 / -J- f (w(z) - S n ) ( j ) 2 (z)f n ^(z)dz. (2.44) 

o r(y) i 

In case (I), it can be readily verified that because / n (oo) = 1 and 

/; < 0, (2.45) 

we have 

/n(0) > f n (x) > / n (00) = 1. (2.461) 



In case (II), we assume w(x) to be not too large so that (2.44) is consistent 
with 

f n {x) > for all x 

and therefore 

/n(0) = 1 > f n (x) > / n (oo) > 0. (2.46II) 

As we shall see, these two boundary conditions (I) and (II) produce sequences 
that have very different behavior. Yet, they also share a number of common 
properties. 

Now, assuming that u'(x) < (section 3) and hence w'{x) < (by (2.18) 
and (2.28)), we have from refs.[3,4] 

Hierarchy Theorem. (I) With the boundary condition f n (oo) = 1, we have 
for all n 

S n+ 1 > £n (2.47) 

I ( w) < at any x > (2 ' 48) 

Thus, the sequences {S n } and {f n {x)} are all monotonic, with 

Si < S 2 < S 3 < ■ ■ ■ (2.49) 

and 

1< h(x) < f 2 (x) < h{x) < ■ ■ ■ (2.50) 

at all finite x. 

(II) With the boundary condition / n (0) = 1, we have for all odd n = 21 + 1 
an ascending sequence 

Ei < S 3 < S 5 < ■ ■ ■ , (2.51) 
but for all even n = 2m, a descending sequence 

S 2 > S A > S 6 > ■■■ . (2.52) 

Furthermore, between any even n = 2m and any odd n = 21 + 1 

82m > 821+1. (2.53) 
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Likewise, at any x, for any even n = 2m 

d_ I hm+\{x) 
dx \ f 2m {x) 

whereas for any odd n = 21 + 1 



< 0, (2.54) 



dx (/ 2 z+iO)) > °' {255) 
The groundstate energy E of the original Hamiltonian 

H = -\i +v(x) (2 - 56) 

is the limit of the sequence {E n } with 

E n = E - S n . (2.57) 
Hence the boundary condition f n (oo) = 1 yields a sequence 

Ei > E 2 > E 3 > ■ ■ ■ > E (2.58) 

with each member E n an upper bound of E, similar to the usual variational 
method. On the other hand, with the boundary condition f n (0) = 1, while the 
sequence of its odd members n = 21 + 1 yields a similar one, like (2.58) with 

E x > E 3 > E 5 > ■■■ > E, (2.59) 

its even members n = 2m satisfy 

E 2 < E 4 < E 6 < ■ ■■ < E. (2.60) 

It is unusual to have an iterative sequence of lower bounds of the eigenvalue 
E. In both cases we have 

lim E n = E. (2.61) 

rwoo v ' 
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3. Proof of u'(x) < 

In this section we shall establish u'(x) < for the case a = 2 [i.e., the one- 
dimensional case of the generalized sombrero potential (1.2)]. The general 
case of arbitrary positive a will be discussed in the Appendix. 
Write the potential V{x) of (1.6) as 

V{x) = g 2 v(x). (3.1). 

When a = 2, 

v( x ) = i(x 2 - l) 2 (x 2 + 2). (3.2) 
The functions So(x) and Si(x) introduced in Section 2 are 

Q 

Sq(x) = j(x 2 - l)v^T2 - - \n(x + v^T2) (3.3) 

and 

Sl (x) = ln(x + 1) + \ ln(z 2 + 2) + i In * + *+>gfl±3 (3.4) 

4 2 2 — x + y3(ar + 2) 

In accordance with (2.17), we find 

, , 3 25x 8 + 150x 6 + 393x 4 + 408rc 2 + 144 

u{x) = 8 FTWTsj (3 - 5) 



with 
and 

Its derivative is 



A = 5x 6 + 10x 4 + 21x 2 + 12 (3.6) 
B = 8V3x(x 2 + l)Vx 2 + 2. (3.7) 



[ } ~ 8 v^+2(« + PVx^T2y [ } 



with 
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Ci = 250rc 17 + 3000x 15 + 16740x 13 + 49880X 11 + 83706x 9 + 77952rc 7 

+34008x 5 + 2880x 3 - 1152a:, (3.9) 

C 2 = 8\/3(322x 12 + 2236x 10 + 6322x 8 + 9672x 6 + 8904rc 4 

+4800x 2 + 1152), (3.10) 

a={x 2 + 2) 2 A (3.11) 

and 

(3 = ( x 2 + 2)*B = 8V3x(x 2 + l)(x 2 + 2) 2 . (3.12) 

Because at all x 

(4800x 2 + 1152)8V3 > 1152xVx 2 + 2, (3.13) 

we have 

CiVx 2 + 2 + C 2 > (3.14) 

and therefore 

u' < 0. (3.15) 

From (3.5), we see that u > and therefore Hierarchy Theorem is applicable. 
The general case when V(x) is given by (1.6) will be discussed in the Appendix. 
As we shall see, for V(x) = g 2 v(x), 

v( x ) = h x 2 - l) 2 ( x 2 + a) (3.16) 

and 

a> a c ^ 0.664, (3.17) 

the corresponding u{x) also satisfies (3.15); i.e. u'{x) < for all x > 0, and 
that ensures the applicability of the Hierarchy Theorem. 
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4. Numerical Results and Discussions 

Throughout this section we denote the n th order iterative solution as 

1p n (x) = lf> (x)f n {x) (4.1) 

with f n (x) the solution of (2.39) and ip (x) related to <f>(x) of (2.10) by 

M*) = (4-2) 

so that 

^o(O) = 1. (4.3) 



4.1 The case a = 2 and g = 1. 

We first discuss the case of a = 2 and g = 1. The resultant series of 
energies are shown in Table 1. The energy series for two different boundary 
conditions converge in different ways. 

I . For the boundary condition / n (oo) = 1 we find that the iterative energy 
sequence 

E > E x > E 2 > E 3 > ■ ■ ■ > E (4.4) 

is given numerically by 

I. 7321 > 1.0163 > 1.0031 > 1.0005 > 1.0001 > 1.0000 > • • • > E = 1. 

(4.5) 

II. For the boundary condition f n (0) = 1, E is still 1.7321, but the iterative 
energy sequence becomes 

Ei > E 3 > ■ ■ ■ > E > ■ ■ ■ > E 4 > E 2 , (4.6) 

and is given numerically by 

1.0163 > 1.0002 >•• >£">•• •> 1.0000 > 0.9981 . (4.7) 

To 4 decimal places £4 is indistinguishable from the exact E = 1. 

In both cases I and II, the zeroth order trial function ^o(x) has a maximum 
near x = 1. However, to the accuracy of Figure 1, for n = 2 ip2(x) is essen- 
tially the same as the exact groundstate wave function ip(x) = e~ x / 4 . The 
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resultant series of wave functions and energies convergence rapidly to the exact 
expressions. It is interesting to notice that although the trial function has its 
two maxima near x = ±1, the iterative series gives the final wave function in 
the shape of the exact groundstate wave function with a single maximum at 
x = 0. The rapid change of the shape of the wave function from the trial one 
to the exact solution shows how the iteration procedure works. 

4.2 The case g = 1. 

Next we discuss the case of g = 1 for different values of a. The results of 
the energy series for the boundary condition II are given in Table 2. 

As proved in Appendix, to ensure the convergence of the iterative series the 
parameter a should satisfy a > a c = 0.664. We have chosen g = 1, a = 3, 2 
and 1.8. The obtained wave functions are shown respectively in Figures 2, 1 
and 3. The behavior of the final function for a < 2 is similar to the one at 
a = 2, namely the resultant wave function has only one maximum at x = 0, 
while the trial function has two maxima near x = ±1. When a = 3 the itera- 
tive wave function retains a similar shape as the trial function with its maxima 
at x near ±1, indicating that the exact groundstate wave function also has its 
maximum at x near ±1, like the trial function when g = 1 and a = 3. 

4.3 The case a = 2. 

Next we discuss the case a = 2 for different values of g. The results of the 
energy series for the boundary condition II are given in Table 3. 

According to (2.12), for a = 2 the allowed values of g is g > 0.866. The 
wave functions for a = 2, g = 3, 1 and 0.88 are shown respectively in Figures 
4, 1 and 5. The behavior of the final function for g < 1 is similar to the 
one at g = 1, namely the resultant wave function has only one maximum at 
x = 0, while the trial function has two maxima near x = ±1. When g = 3 
the iterative wave function retains a similar shape as the trial function with 
its maxima at x near ±1, indicating that the exact groundstate wave function 
also has its maximum at x near ±1, like the trial function when g = 3 and 
a = 2. 
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Table 1. Eigenvalues of groundstates for g = 1 and a = 2 







Ei 


E2 


^3 


E A 


E 5 


I 


1.7321 


1.0163 


1.0031 


1.0005 


1.0001 


1.0000 


II 


1.7321 


1.0163 


0.9981 


1.0002 


1.0000 


1.0000 



Table 2. Eigenvalues of groundstates for g = 1, 
using II with / n (0) = 1 as the boundary condition 



a 


£0 


£1 


E2 


£3 


£4 


£5 


1.8 


1.6733 


0.9558 


0.9418 


0.9432 


0.9431 


0.9431 


2 


1.7321 


1.0163 


0.9981 


1.0002 


1.0000 


1.0000 


3 


2.0000 


1.2974 


1.2602 


1.2659 


1.2651 


1.2652 



Table 3. Eigenvalues of groundstates for a = 2 
using II with / n (0) = 1 as the boundary condition 
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E 


Ei 


E 2 


E 3 


E A 


E 5 


0.88 


1.5242 


0.8633 


0.8517 


0.8528 


0.8527 


0.8527 


1 


1.7321 


1.0163 


0.9981 


1.0002 


1.0000 


1.0000 


2 


3.4641 


2.6934 


2.6375 


2.6465 


2.6455 


2.6456 


3 


5.1962 


4.5786 


4.5562 


4.5591 


4.5589 


4.5589 
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Appendix 

In this Appendix we give the proof of the convergence of the iterative series 
for the generalized double-well potential, eq.(1.6) 

V(x) = g 2 v(x) and v(x) = hx 2 - l) 2 (x 2 + a), (Al) 

where a is an arbitrary constant. 

Following the steps described in Section 2, from (2.4), (2.7) and (2.17) we 
have 

/ \ a-8y / x 2 + a P 

U[X) ~ 8(x 2 -l) 2 (x 2 + a) 2 [ ] 

where 

a = 15x 6 + 6(3a - l)x 4 + (8a 2 + 12a + l)x 2 + 8a 2 + 2a (A3) 

and 

= VI + a x(3x 2 + 2a- 1). (A4) 
Integrating (2.4) and (2.7) gives 

1 , i—z x q /a 1 



So(x) = -x(\/x 2 + a) 3 — ( - + - jxVx 2 + a 

~ a (f + \) ln ( x + V^ 2 + a) ( A5 ) 



and 



= ln(* + + + \ in f±Ifg + g + g (A6) 

^ Va + lv £ + a + a — x 

The corresponding trial function (2.10) satisfies the Schroedinger equation 
(2.16) with w{x) expressed by (2.18) and (2.19). 

According to the Hierarchy Theorem proved in Ref.[3,4] the iterative series 
is convergent if the potential w{x), defined by (2.18), satisfies the conditions 
w{x) > and w'(x) < 0. As we shall prove, when 

a > a c ^ 0.664, (A 7) 
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the potential u{x), defined by (2.17), satisfies u{x) > and u'(x) < 0. 
Furthermore, when g satisfies (2.12), i.e. 



a > a 9 = TTo (^-8) 



we have T > from (2.13); therefore g{x) > and g'{x) < 0. Combining 
these results we have 

w(x) > and w'(x) < when a > a m = max(a c , a g ). (^4-9) 

This condition ensures the convergent iterative solution for the groundstate of 
the generalized double-well potential (A.l) (i.e. (1.6)), in accordance with the 
Hierarchy theorem. 

Range of parameter a 

As mentioned above, to obtain a convergent iterative series the potential 
u{x) in (A. 2) should satisfy 



u(x) > and u(x) < 0. 

For a > 0, introducing 

7± = a ± 8 V x 2 + a (3 
with a and (3 given by (A.3)-(A.4), we have 

7+7 _ ee (x 2 - 1) 2 7 = a 2 - 64(x 2 + a) (3 2 . 

Thus, u{x) can be expressed as 

7- 

ii = 

8(x 2 + a) 2 (x 2 - l) 2 

and 

u = I m 

8(x 2 + a) 2 V 7+ A 

Write 

7 = 7(a, x) = ^i(x)(15x 4 + 36ax 2 ) + ^2 

where 

gi(x) = I5x 4 +I8x 2 - 1 
16 



A10) 
All) 

A12) 

A13) 

A14) 

A15) 
A16) 



and 

g 2 = g 2 (a, x 2 ) = 4(141x 4 + 90x 2 + l)a 2 + 32(9x 2 + l)a 3 + 64a 4 > 0. (A17) 

Consider only the first quadrant in (x, a)-plane with x, a positive. Define 
xq to be the point satisfying gi(xo) = 0. For x > xq, g\{x) > and therefore, 

7 > 0. For x < xq, gi(x) < and 



d_ 

da W 



(^) > for a > 0. (A18) 



Since (7/a 2 ) — > —00 when a — > and (7/a 2 ) — > 00 when a — > 00, at a fixed 
x < xo by varying the parameter a from to 00, 7(0, x) = only once. We 
define a curve (G) in x, satisfying 7(0, x) = and a curve (5) satisfying 
/?(a, x) = as shown in Fig. Al. It can be seen that above the curves (B) 
and (G), f3 and 7 are both positive. The curve (B) is above the curve (G). 
Therefore above the curves (B) 7+, 7 and u are all positive. Below the curve 
(B) (3 is negative; therefore 7 and u are positive. Thus we can conclude that 
u{x) > everywhere for a > 0. 
By differentiating u, we have 

(ai)a = U ' {X) = 8(x 2 + a)3(x 2 -l) 3 ' (j4 ' 19) 



where 
with 



7_ = a - 8Vx 2 + a (3 (A20) 



a = [x 4 + (a — l)x 2 — aja' — [8x 3 + 4(a — l)x]o; 
= (-30rc 9 - 6x 7 - 42x 5 + Ux 3 ) + (-42x 7 - 162x 5 + 18x 3 - 6x)a 

+ (-48x 5 - 144x 3 )a 2 + (-16x 3 - 48rc)a 3 (A21) 

and 

= ( x 2 + a){x 2 - 1)0' - [7x 3 + (4a - 3)x]P 
= \/aTT[(-12x 6 +6x 4 -2x 2 ) + (-15x 4 -2x 2 + l)a+(-6x 2 -2)a 2 |. (A22) 
Defining 



7± = a ± 8Vx 2 + a /3 (A23) 
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we have 

7+7- = O 2 - 1) 3 7 

= a 2 - 64(> 2 + a) (3 2 . (A.24) 

Therefore 

m'= c , * (A.25) 

8(x 2 + a) 3 v 7+ / v 7 

It is straightforward to show after some derivation that 7 is a polynomial of x 2 
and a. Defining 

7 = £ f A x 2A (A26) 

A=0 

where f A are polynomials of a: 

f = a 3 [64- 192a + 256a 3 ] 

f 1 = a 2 [-228 - 1152a 2 + 1536a 3 ] 

f 2 = a[168 + 1068a- 960a 2 + 3648a 3 ] 

f 3 = 60 - 504a + 4500a 2 + 4992a 3 

f 4 = -180 + 8568a + 4644a 2 

f 5 = 3060 + 2520a 

f 6 = 900. (A27) 

For small a, it is convenient to introduce 

z = x 2 /a, i.e. x 2 = az. (A. 28) 

We can plot three curves on the (z, a)-plane : the curve (A)\ a = 0, (B)\ 
(3 = and (G): 7 = 0. The regions above the three curves correspond to 
a < (3 < and 7 > 0, respectively (see Fig. A2 for details). Above the 
critical point C, i.e. a > a c = 0.664 we have u' < from (A. 23) and (A. 25). 
Therefore this is the condition for u > and u' < 0. 
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Figure Caption 

Fig. 1 Trial Function ipo(x) and Groundstate Wave Function ip{x) 
for g = 1 and a = 2. 

Fig. 2 Trial Function ^o(a^) and Groundstate Wave Function ip{x) 
for g = 1 and a = 3. 

Fig. 3 Trial Function ^o(^) ar| d Groundstate Wave Function ip{x) 
for g = 1 and a = 1.8. 

Fig. 4 Trial Function ^o(a^) and Groundstate Wave Function ip(x) 
for g = 3 and a = 2. 

Fig. 5 Trial Function ^o(^) ar| d Groundstate Wave Function ip{x) 
for g = 0.88 and a = 2. 

Fig. Al Curves for /? = and 7 = 0. 

Fig. A2 Curves for a = 0, (3 = and 7 = 0. 
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Fig. 3 Trial Function ipo( x ) ar| d Groundstate Wave Function ip(x) 

for g = 1 and a = 1.8. 
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Fig. 4 Trial Function ^o(x) and Groundstate Wave Function ip(x) 

for g = 3 and a = 2. 
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Fig. Al Curves for (3 = and 7 = 0.. 
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